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Abstract 

Owing to the structure of the Gaussian multiple-input multiple-output (MIMO) broadcast channel (BC), 
associated optimization problems such as capacity region computation and beamforming optimization are 
typically non-convex, and cannot be solved directly. One feasible approach to these problems is to transform 
them into their dual multiple access channel (MAC) problems, which are easier to deal with due to their 
convexity properties. The conventional BC-MAC duality is established via BC-MAC signal transformation, and 
has been successfully applied to solve beamforming optimization, signal-to-interference-plus-noise ratio (SINR) 
balancing, and capacity region computation. However, this conventional duality approach is applicable only to 
the case, in which the base station (BS) of the BC is subject to a single sum power constraint. An alternative 
approach is minimax duality, established by Yu in the framework of Lagrange duality, which can be applied to 
solve the per-antenna power constraint problem. This paper extends the conventional BC-MAC duality to the 
general linear constraint case, and thereby establishes a general BC-MAC duality. This new duality is applied to 
solve the capacity computation and beamforming optimization for the MIMO and multiple-input single-output 
(MISO) BC, respectively, with multiple linear constraints. Moreover, the relationship between this new general 
BC-MAC duality and minimax duality is also presented. It is shown that the general BC-MAC duality offers 
more flexibility in solving BC optimization problems relative to minimax duality. Numerical results are provided 
to illustrate the effectiveness of the proposed algorithms. 



I. Introduction 

In a Gaussian multiple input multiple output (MIMO) broadcast channel (BC), the base station (BS) 
equipped with multiple transmit antennas sends independent information to each of multiple remote 
users, which are equipped with multiple receive antennas. In the past decade, a great deal of research 
has been focused on the characterization of optimal transmission schemes for the MIMO BC [l]-[6]. 
Due to the coupled structure of the transmitted signals, the optimization problems associated with the 
BC are usually non-convex. The key technique used to overcome this difficulty is to transform the 
BC problem into a convex multiple access channel (MAC) problem via a so-called BC-MAC duality 
relationship. Up to now, two types of BC-MAC duality have been proposed as follows: 

1) Conventional BC-MAC DuaUty ( [1], [3], [7], [8]): 

Under a single sum power constraint, the capacity region (or signal-to-interference-plus-noise- 
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ratio (SINR) region) of a BC is identical to that of a dual MAC under the same sum power 
constraint. The channel matrix associated with the dual MAC is the conjugate transposed channel 
matrix of the BC, and the noise co variance matrices of both channels are identity matrices [1], 
[9]. 

2) Minimax Duality ( [5], [10], [11]): 

The sum rate maximization problem of a BC with multiple linear constraints has the same 
solution as the dual MAC minimax optimization problem. The channel matrix of the dual MAC 
is the conjugate transposed channel matrix of the BC, and the noise covariance matrix of the dual 
MAC is an unknown variable of the minimax optimization problem [10]. 
The conventional BC-MAC duality was first observed by Rashid-Rarrokhi et al. [1], and applied to 
solve the sum power minimization problem for a BC with SINR constraints. Several different methods 
have been developed independently to prove the conventional BC-MAC duality. The proof in [1] is 
based on the equivalent transformation that maps the SINR of the MAC to that of the BC. Vishwanath 
et al. [12] proved the conventional BC-MAC duality by presenting the explicit transformation between 
the transmit covariance matrix of the BC and that of the MAC, and applied this duality to solve the 
sum capacity problem. Both proofs here are based on the reciprocity relation [13] between the BC 
and its dual MAC. Another proof based on the Karash-Kuhn-Tucker (KKT) conditions was given by 
Visotsky and Madhow [14]. The conventional BC-MAC duality has been widely applied to many BC 
problems. Schubert and Boche [15] solved an BC SINR balancing problem, which is to maximize the 
minimal SINR among all the users under a sum power constraint, via transforming the problem into 
its dual MAC problem. The conventional BC-MAC duality was also used to show that a dirty paper 
coding (DPC) [16] is a sum-capacity achieving strategy by Viswanath and Tse [8]. Moreover, the entire 
capacity region for the MIMO BC channel can be obtained via the conventional BC-MAC duality [3], 
[4], [17]. However, the conventional BC-MAC duality is applicable only to the case in which the BS 
of the BC is subject to a single sum power constraint. 

On the other hand, the sum-capacity for the MIMO-BC was also studied by Yu and Cioffi [11] via 
minimax optimization. The minimax duality was proposed by Yu [10], where the conventional BC- 
MAC duality and the minimax duality are unified in the framework of Lagrange duality. However, only 
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the sum capacity is considered in [10]. Furthermore, Yu and Lan [5] extended the minimax duality to 
solve the capacity region computation problem and beamforming problem for the BC with per-antenna 
power constraint. The proofs of the minimax duality in [5] and [10] are based on Lagrange duality. 
Compared with the conventional BC-MAC duality, the minimax duality can be applied to the case 
with multiple linear constraints. However, since the dual MAC problem has a minimax form, and the 
noise covariance matrix of the dual MAC is unknown, high-efficiency algorithms, such as the iterative 
water- filling algorithm [18], cannot be applied. 

A. Overview of the Main Results 

The purpose of this paper is to establish the general BC-MAC duality via the BC-MAC SINR 
transformation, and unify the BC-MAC duality and the minimax duality in the framework of the 
reciprocity relationship between the BC and the MAC. By introducing several auxiliary variables and 
applying the general BC-MAC duality, the primal BC problem with multiple transmit covariance 
constraints is transformed into its dual MAC problem with a single sum power constraint and can 
be efficiently solved via the existing algorithm for its dual MAC as the MAC problem has a convex 
structure that is easier to handle. 

In this paper, we first consider a MIMO BC with a single general linear constraint. Relying on the 
BC-MAC transformation, we prove that the capacity region of the BC is the same as that of its dual 
MAC with a single weighted sum power constraint, which we term the general BC-MAC duality. The 
channel matrix of the dual MAC is the transposed channel matrix of the primal BC, and its noise 
covariance matrix is the coefficient of the linear constraint instead of being an identity matrix as in the 
conventional BC-MAC duality. 

To exploit the general BC-MAC duality, the weighted sum rate maximization problem for the BC 
with multiple linear constraints is transformed into a single linear constraint problem by introducing 
several auxiliary variables. Though the rate maximization problem for the BC is a non-convex problem, 
we show that the KKT conditions are sufficient for optimality, and show that the subgradient-based 
algorithm converges to the optimal solution. A method for obtaining the subgradient is also given. 
The relationship between the general BC-MAC duality based solution and the minimax duality based 
solution [5] is explored. We show that the new method to handle multiple-constraint is equivalent to 

September 24, 2008 DRAFT 



4 



that of the minimax duality. But since the general BC-MAC duality based method solves the multiple 
constraint optimization problem in a decoupled manner, the new result has more flexibility to apply 
the existing algorithms for the MAC, while the minimax duality does not. Moreover, we discuss the 
weighted sum rate maximization problem with a convex but nonlinear constraint, and develop a new 
iterative algorithm to solve this optimization problem. In addition to the weighted sum rate maximization 
problem, the proposed method is also applied to solve the beamforming problem in a multiple-input 
single-output (MISO) BC with multiple linear constraints. 

B. Organization and Notation 

The remainder of this paper is organized as follows. The system model is described in Section II. 
Section III presents the general BC-MAC duality, where the transmit covariance matrix of the BC 
is subject to a linear constraint. The capacity region computation problem of the BC with multiple 
constraints or a single nonlinear constraint is studied in Section IV. The method to cope with multiple 
linear constraints is also applied to solve the beamforming problems in Section V. Several numerical 
results are provided in Section VI to illustrate the effectiveness of the proposed methods. Finally, 
Section VII concludes the paper. 

Throughout this paper, we use boldface upper and boldface lower case letters for matrices and vectors, 
respectively. (■)^ denotes the matrix conjugate transpose operation, and tr(-) denotes the matrix trace 
operation. E(-) denotes the expectation operation for random variables. / denotes an identity matrix. 

II. System Model 

We consider a MIMO BC system shown in Fig. 1 (a), where the BS intends to send independent 
information streams to each of K remote users. The BS has Nt transmit antennas and each user has 
Nr receive antennas. The signal received by the zth user is modeled as follows: 

= HiX + Zi, i = !,■■■ ,K (1) 

where the A^f x 1 vector x denotes the transmit signal at the BS, Hi denotes the Nr x Nt channel matrix 
from the BS to the ith user, denotes the receive signal at the ith user, and Zi is the noise vector. 
The entries of Zj are modeled as independent and identically distributed complex Gaussian random 
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variables with mean zero and variance af. The transmit signal covariance matrix of the BS is defined 
as Q := 'E(xx^). In this paper, we assume that the channel knowledge is perfectly known at both the 
BS and the users, i.e., is perfectly known at the transmitter and the receivers. 
Since the transmit signals for different users are independent, we have 



K 

X 

i=l 

where Xi denotes the transmit signal intended for the ith user. Furthermore, to fully utilize the spatial 



Y^Xi (2) 



diversity of the multi-antenna system, a spatial multiplexing scheme is also applied, which means that 
the data intended for each user is further divided into N substreams, where = mm{Nt., N^) [13]. 
Thus, the transmit signal for the iih user can be expressed as follows: 

N 

Xi ^ ^ ^ij^ij 

where bij is a complex scalar variable with pij := E(|6jjp), representing the information signal of 
the jth data substream of the ith user, and Uij denotes the corresponding beamforming vector with 
\\uij\ \ = 1. Combining (2), (3) and the definition of the transmit covariance matrix Q, we have 

K N 

i=i j=i 

A. Nonlinear Encoding and Decoding Strategy 

It has been shown that the DPC scheme is a capacity achieving scheme for the Gaussian MIMO 
BC [17]. With the DPC scheme, the information for different users is encoded in a sequential manner. 
Without loss of generality, we assume that the encoding order is identical to the index order, i.e., the 
data substream 6i i is encoded first, 6i 2 is next encoded, and so on. According to the DPC scheme, 
the latter encoded data stream has non-causal information about its former encoded data streams, and 
thus the interference caused by the former data streams' transmission can be completely removed by 
the DPC scheme. Thus, the rate achieved by the ith user can be expressed as 

where Qj := E{xixf'). 
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(a) BC, Zi ~ Af{0, all), QA<P 
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(b) Dual MAC, z ~ Af{0, A), J2f=i ^MQT^) < P 



Fig. 1. The system models for the primal MIMO BC and the dual MAC. 



At the receiver side, for each user, successive interference cancellation (SIC) and the linear minimum 
mean square error (MMSE) filter are adopted to decode the corresponding information. With SIC, the 
first data stream is decoded by treating all the other streams as interference; then the signal from the 
first data stream is subtracted from the received signal, and the second data stream is decoded next, 
and so on. Thus, the mutual information between the BS and the ith user can be expressed as 



I{hiX,yi\h 



i.li 



1 bi,N-i 



)• (6) 



Moreover, since the MMSE receiver is information-lossless [19], each term in (6) is achievable with 
the MMSE receiver. Thus, the MMSE-SIC receiver can achieve the capacity of the MIMO system. The 
receive beamforming vector for the jth data substream at the zth user is denoted by the A^^ x 1 vector 
Vij. Thus, the SINR of the jth data substream at the ith user receiver can be written as 



SINK. 



Since the achievable rate of each data substream depends on its SINR, we have 



(V) 



N 



n = ^log(l + SINR 



(8) 
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B. Linear Encoding and Decoding Strategy 

Although the nonlinear strategies DPC and SIC are capacity-achieving schemes, they are difficult to 
implement in practice. A straightforward scheme for transmission is beamforming without DPC at the 
transmitter side and SIC at the receiver side. In the linear strategy, the transmit and receive beamforming 
vectors for the jth data substream of the ith user are still denoted by Uij and Vij, respectively. Thus, 
the corresponding SINR can be expressed as 

CTNrW - PijKl^i^iJ^ (Q. 

Lfc=l L/=l {k,l)^{^d)Pk,l\W,jHiUk,l\^ + af 

Since the method developed in the present paper is applicable to both linear and nonlinear strategies, 
we mainly focus on the nonlinear strategy, and adopt the SINR definition (7) throughout the paper. 

C. General Linear Transmit Covariance Constraint 

In the aforementioned literature, the transmit covariance matrix is subject only to a sum power 
constraint or/and a per-antenna power constraint. In this paper, we consider a general linear transmit 
covariance constraint as follows: 

tr(QA) < P (10) 

where A is an Nt x Nt matrix, and P is a constant. If the matrix A is chosen to be an identity 
matrix, then the constraint (10) is reduced to the sum power constraint; if the matrix A is chosen to 
be the diagonal matrix having all diagonal elements being zero except the jth element being 1, then 
the constraint (10) is reduced to the jth antenna power constraint. In cognitive radio networks, we 
choose A = hh^ , where h is the channel response from the secondary transmitter to the primary 
receiver, and the constraint is reduced to the interference power constraint for protecting the primary 
users [20]-[23]. Therefore, the constraint (10) can be viewed as a generalized linear constraint. 

D. Objective Functions 

In this paper, we consider several scenarios with different objectives: MEMO BC capacity region 
computation, SINR balancing, and power balancing. 
1 ) MIMO BC Capacity Region Computation 
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Any boundary rate tuple of the Gaussian MIMO BC capacity region can be obtained by solving the 
weighted sum rate maximization problem with some given user rate weights. Therefore, the capacity 
region computation problem is formulated as follows: 

K 

max > WiTi (11) 

where Wi is the positive weight of the zth user, and Q > denotes the semidefiniteness constraint. By 
varying the values of the weight WjS, the entire capacity region of the MIMO BC can be obtained. 

2) SINK Balancing 

The aim of the SINR balancing problem is to maximize the minimal ratio between the achieved 
SINR and the target SINR among all the data substreams. Mathematically, the optimization problem 
is formulated as 

• SINRi,- 

max mm (12) 

<3>o Vij 'jij 

where jij is the target SINR for the jth data substream of the ith user. Conventionally, the SINR 
balancing problem considers the MISO case [15], i.e., Nt = 1. It has been shown in [15] that the ratios 
of all the data substreams are equal, when the optimal solution is achieved. Therefore, the problem is 
termed the SINR "balancing" problem. Note that (12) is equivalent to the following form 

max a 

subject to SINRjj > a'jij, V 
In the sequel, we will use (13) as the optimization problem for the SINR balancing problem instead 
of (12), since it is easier to write the Lagrange function of (13). 

3) Power Balancing 

In this case, the system has several different power requirements, and the objective is to minimize the 
maximal ratio between the transmit power and the power requirement. Mathematically, the optimization 
problem can be expressed as 

tr(QA,) 

mm max (14) 

Q>0 Vi Pi 
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where Pi is the ith power requirement, and Ai is an Nt x Nt matrix for the ith power requirement. 
Similarly to the SINR balancing function, the problem (14) can be transformed into 

min a 

Q>0,a ^^^^ 

subject to tr(QAj) < aPi, Mi. 
If there is a single power requirement, and the corresponding matrix A is an identity matrix, then (14) 
reduces to the power minimization problem. 

In this paper, we will consider these optimization problems with several general transmit covariance 
constraints or SINR constraints. 

III. General BC-MAC Duality 

In this section, we establish the general BC-MAC duality under a single linear transmit covariance 
constraint. We start with the SINR balancing problem expressed as follows: 

max a 

Q>0,a 

subject to SINRij > a-fij, \f i,j (1^) 
tr(QA) < P. 

The MIMO BC SINR balancing problem (16) is a non-convex optimization problem due to the non- 
convex SINR constraints. Although it has been shown in [5], [24] that the SINR constraint under the 
MISO scenario can be transformed into the second order cone (SOC) form, the transformation is not 
applicable to the MIMO case due to the essentially non-convex property of the MIMO SINR constraints. 
Hence, the problem (16) is still an open problem, and cannot be solved via existing methods. However, 
we can establish a new MAC called the dual MAC, and formulate a dual MAC problem of the primal 
problem (16) such that it shares the same solution as its primal problem (16). 

Definition 1: The dual MAC of the primal BC in (1) has the conjugate transposed channel matrix 
of the BC, i.e., the channel matrix of the dual MAC from the ith user to the BS is , and the noise 
covariance matrix at the BS is the matrix A instead of the identity matrix, which is shown in the Fig. 
1 (b). 
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The corresponding dual MAC optimization problem is expressed as follows: 

max a. 

subject to SINrJJ^ > 07^^, V z,j (17) 

i=\ 

where Q'f^^ denotes the transmit signal covariance matrix of the zth user, SINrJ"^ denotes the SINK 
of the jth data substream of the ith user, and the superscript '(m)' indicates that the corresponding 
variables are for the dual MAC. In this dual MAC, the MMSE-SIC scheme is applied, which means 
that the data streams of the dual MAC are decoded in a sequential manner. In the dual MAC problem 
(17), the decoding order at the BS is converse to the encoding order of the DPC at the primal BC. 
Accordingly, the SINR^™'' can be expressed as 

where g^j is the power allocated to this data substream, and itj ^ and Vi j denote the corresponding 
receive and transmit beamforming vectors, respectively. While it may be somewhat confusing at first 
that the beamforming vectors of the dual MAC share the same notation with the beamforming vectors 
of the primal BC in (7), it will become clear in the following that the optimal beamforming vectors of 
the primal BC are identical to those of the dual MAC. 

The relationship between the primal problem (16) and the dual MAC problem (17) is summarized 
in the following proposition. 

Proposition 1: The optimal solutions of the primal problem (16) and the dual MAC problem (17) 
are the same. 

Proof: First, we will prove by contradiction that the SINR constraints for problem (16) hold 
with equality when the optimal solutions are achieved. If SINRj > a*7jj with a* being the optimal 
solution of problem (16), then we can reduce part of the power pi j and distribute it to all the other data 
substreams, thereby increasing the objective value a without violating the constraints. This contradicts 
a* being the optimal solution. A similar argument holds for the problem (17). 
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We next prove that if a is achievable for the problem (16), then it is also achievable for the problem 
(17). Assume that when a is achieved, Uij and ■Uj j are the corresponding beamforming vectors for 
transmitter and receiver, respectively, and pij is the power allocated to the jth data substream of the 
ith user. For the dual MAC problem (17), we can choose Vij to be the transmit beamforming vector 
for the user side, and Uij to be the receive beamforming vector at the BS. The power allocated to the 
jth data substream of the ith user of the dual MAC is assumed to be qij, which can be obtained by 
setting SINrJ™' = SINR^j = a-fij, i.e.. 





Pi,j\ 


\vfjHiUij\ 


2 




\v^jHiUkj\ 


? + El 


j+1 Pi,l 


\v^jHiU^\ 





( Er=\ El, -quM^Kivl^H, + ECi %iHfv.,iv^,H, + a) u 



(19) 



By rearranging (19), we can list all the equations related to the SINR as follows: 

K N N 

k=2 1=1 1=2 

K N N 

Pi,2{ul^{qi,iH^Vi^ivliHi + A)ui^^ = qi^2[^J2^PkAvf2HiUk,i\'^ +^Pi,i|i;f2-H"i'Ui,zP + f^i 

k=2 1=1 1=3 
K-1 N N-1 

Pk,n(u^^n[Y1 Y1 ^k,iH^Vk,iv^^iHk + qK,iH^VK,iVK,iHK + A^Uk,n^ = qK,NcrK- 

k=l 1=1 1=1 

By adding the above equations together, we have 

K N K N 

i=l j=l i=l j=l 

From the power constraint tr(QA) < P of the problem (16), we have Eli Ej'=iPi,j'^i^j^'^iJ — P- 
Combining this with (20), we have Eli E!j=i ^ which means that the weighted sum power 
constraint Eli (^i^^iQi^^) < -P of the problem (17) is satisfied. Hence, a is achievable for the problem 
(17). 

Contrarily, we can prove that if a is achievable for the problem (17), then a is also achievable for 
the problem (16). 

The proof of Proposition 1 follows. ■ 
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If we assume that the optimal solution for both problems (16) and (17) is a*, then it means that the 
point {SINRjj = Q;*7jj, V i,j} is achievable for the primal BC and the dual MAC with corresponding 
constraints, respectively. Therefore, under the general linear constraint tr(QA) < P, the primal BC can 
achieve the same SINR region as the dual MAC, which is subject to a weighted sum power constraint. 
From an information theoretic perspective, according to (8), the rate point {r = [ri,-- - ,rx]} is 
achievable if the SINR point {SINRi Vz, j} is achievable under corresponding constraints. Thus, we 
have the following corollary. A rigorous proof is provided in Appendix A. 

Corollary 1: The capacity region of the primal BC under the constraint tr(QA) < P, is equal to the 
capacity region of its dual MAC with a single weighted sum power constraint Xl^i (^i^^iQi^^) ^ P- 

Remark 1: By setting A in (16) to be an identity matrix, and assuming that af = 1 for all users, the 
general linear power constraint becomes a sum power constraint, and the noise covariance at the BS of 
the dual MAC is reduced to an identity matrix. This is precisely the same as the conventional BC-MAC 
duality. Thus, the new BC-MAC duality can be viewed as a generalization of the conventional BC-MAC 
duality. The proof of the duality is based on the special BC-MAC reciprocity relationship, instead of 
the Lagrange duality used in [5], [10]. Note that since the SINR constraints in (16) is not convex, 
the Lagrange duality gap between (16) and (17) may not be zero. Therefore, Proposition 1 cannot 
be proved through the use of Lagrange duality. But for the capacity region problem, the objective 
function is concave in the signal covariance matrices and convex in the noise covariance matrices. 
Hence, the Lagrange duality gap is zero and the Lagrange duality can be applied for the proof of the 
BC-MAC capacity duality. From this perspective, the reciprocity relationship is more fundamental than 
the Lagrange duality for the BC-MAC duality. 

Remark 2: For a given set of transmit covariance matrices of the dual MAC, we can obtain the 
corresponding transmit covariance matrix of the primal BC to achieve the same value of a by using 
the method giving in the proof of Proposition 1 . The detailed MAC-BC covariance matrix transformation 
algorithm is provided in Table I. Similarly, the BC-MAC covariance matrix transformation algorithm 
is readily obtained. Furthermore, the proof of Corollary 1 presents a MAC-BC covariance matrix 
transformation such that the primal BC and its dual MAC achieve the same rate tuple. The detailed 
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TABLE I 

The MAC-BC covariance transformation for SINR equivalence. 
MAC-BC covariance matrix transformation I 

1. Apply eigenvalue decomposition to Q^™' — ViAiVi with Vij being the jth column of Vi 

and Qij being the jth diagonal element of Ai, 

2. For each data substream, apply the MMSE algorithm to compute the receive beamforming vector Ui,j, i.e., 

3. According to (19), compute V 

TABLE II 

The MAC-BC covariance transformation algorithm for capacity equivalence. 
MAC-BC covariance matrix transformation II 

1. Define ff^ = HkA~^/\ Mt, = {l^ Efl.+ i fff Q^'"^^,) and 

b, = (i + h,{y.U+iQj)h") 

2. fov k = ,K 

3. calculate the SVD decomposition of M^^''^ H kB^^'^ = RkDkLk 
5. end for 



algorithm is presented in Table II. 

Moreover, although Proposition 1 is for the nonlinear scheme, it is also applicable to the linear 
transmission scheme. 

Corollary 2: Under the linear transmit strategy, the achievable SINR region of the primal BC under 
the constraint tr(QA) < P, is equal to the achievable SINR region of its dual MAC with a single 
weighted sum power constraint ^f^i aftr{Q^[^^) < P. 

In the following, we will show how to use this general BC-MAC duality to solve various BC 
optimization problems with multiple transmit covariance constraints. 
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IV. Capacity Computation 

In the preceding section, the duality between the MIMO BC and the dual MAC with a general linear 
constraint was presented. In this section, by exploiting this duality, we will compute the capacity region 
of the MEMO BC. According to the discussion in Section II-D, the capacity of the BC can be obtained 
by solving the weighted sum rate maximization problem. For simplicity, we assume in the sequel that 
af = 1 for all the users. 

A. Single Linear Transmit Covariance Constraint 

We first consider the weighted sum rate maximization problem for the BC with a single linear 
constraint, which is formulated as follows: 

K 

max y WiVi 

tr (21) 

subject to tr(QA) < P 

where A is a constant matrix, and P is a constant. The problem (21) is a non-convex problem, and 
thus cannot be solved directly. According to Corollary 1, the problem (21) is equivalent to its dual 
MAC problem as follows: 

K 

E(m) 
WiTl ' 

Wr'">0 i-1 

; (22) 



subject to ^tr(QS™^) < P 



where r)" := loe;^ — — — - denotes the achievable rate of the ith user. By solving the 
problem (22) via the interior point algorithm [25], the optimal solution for the problem (21) can be 
obtained via a MAC-BC transmit covariance matrix transformation algorithm. 

In the following, we present an important property of the problem (21), which will be used in the 
case with multiple transmit covariance constraints. We first list the KKT conditions of the problem 
(21) as follows: 



dQi (23) 
A(^tr(gA) -P) = 
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where A is the Lagrange multiplier, and is the Lagrange multiplier associated with the constraint 
Qj > 0. In general, the KKT conditions are only necessary for a solution to be optimal for a non- 
convex problem. However, for the problem (21), it is shown in the following proposition that the KKT 
conditions are also sufficient for optimality. 



Proposition 2: The KKT conditions (23) are sufficient for a solution to be optimal for the problem 
(21). 



Proof: According to the Corollary 1, the problem (21) is equivalent to its dual MAC problem 
(23). We now assume that Q satisfies the KKT conditions in (23) and achieves the weighted sum rate 
R. Then, by the BC-MAC transmit covariance matrix transformation, we can obtain a set of q!'"^^ for 
the problem (22) to achieve the same R. We next assume that Q is an optimal solution of the problem 
(21) with the optimal weighted sum rate R, where R > R. Thus, we can obtain the optimal solution 
of the problem (22) Qj^™"* by MAC-BC transmit covariance matrix transformation. 

It is well known that the objective function of (22) is a convex function. Hence, we have Q* := 
O!™^ +^ (^Qi™'' ~ Qi™^ j ' where < t < 1, is a better solution than Q,-™'' for the problem (22). Through 
the MAC-BC transmit covariance matrix transformation, we transform the dual MAC solution Q* into 
its corresponding BC solution Q*. Since the MAC-BC transmit covariance matrix transformation is 
continuous, we can always find a t such that ||Q — Q*|| < e for a given e > 0. Therefore, Q is not the 
local optimal solution, which is contradicted with the KKT condition di^f^-^WiTi) / dQi = XA + ^'j. 
The proof thus follows. ■ 

B. Multiple Linear Transmit Covariance Constraints 

We now consider the weighted sum rate maximization problem with multiple linear constraints as 



follows: 



K 



max 

Q>0 




i=l 



subject to tr(QAi) < Pi 



(24) 



tr(QA2) < P2 
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where Ai, i = 1,2, is an Nt x Nt constant matrix, and Pj, z = 1,2, is a constant. For convenience 
of description, we discuss only the case of two linear constraints, though our method can be easily 
extended to the case of an arbitrary number of linear constraints. 

Since the objective function is non-concave in Q, the problem (24) is not convex, and thus cannot 
be solved directly. In [3], [4], the sum capacity of the MIMO BC with a single sum-power constraint 
was studied. Based on the conventional BC-MAC duality [12], the BC problem was transformed into 
its dual convex MAC problem with a single sum-power constraint. However, the problem (24) is with 
multiple constraints, and thus the conventional BC-MAC duality cannot be applied. To exploit the 
general BC-MAC duality, we can transform the problem (24) into the following problem with a single 
constraint: 



where Ai and A2 can be viewed as auxiliary variables. The relationship between the problem (24) and 
the problem (25) can be summarized as follows. 

Proposition 3: The optimal solution of the problem (25) is an upper bound on that of the problem 



Proof: If Q is feasible for the problem (24), then it is also feasible for the problem (25). Therefore, 
the feasible region of the problem (24) is a subset of that of the problem (25). The proof follows 



Furthermore, we can prove that the upper bound is tight, i.e., the optimal solution of the problem 
(24) achieves the upper bound g{\i, X2) for some Ai and A2. Thus, we have 

Proposition 4: The optimal solution of the problem (24) is equal to that of the problem minAi,A25'('^i) ■^2)- 
Proof: The KKT condition of the problem (24) can be listed as follows: 



K 




(25) 



subject to Aitr(QAi) + A2tr(QA2) < AiPi + A2P2 



(24). 



immediately. 



d J2i=i Win 

dQi 



/iiAi + /i2A2 + fij, Vi 




(26) 
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where fii and fj,2 are the Lagrange multipliers with respect to the two constraints, respectively, and i7j 
is the Lagrange multiplier associated with the constraint Qi > 0. When the optimal solution of the 
problem (24) is achieved, we assume that the corresponding optimal variables are Q*, /i*, /i2, and Ct*. 
We now list the KKT conditions of the problem (25) as follows: 

OQi (27) 
A(Aitr(gAi) + A2tr(QA2) - fxiPi - /i2P2) = 
where A is the Lagrange multiplier, and Tj is the Lagrange multiplier associated with the constraint 
> 0. If we choose Q = Q*, A = 1, Ai = fil, A2 = fJ^l, and Tj = Q*, the KKT conditions (27) 
are satisfied. Since the problem (25) is the weighted sum rate maximization problem with a single 
linear constraint, according to Proposition 2, the solution is the optimal solution of the problem (25). 
Combining this with Proposition 3, the proof follows. ■ 
According to the general BC-MAC duality discussed in Section III, the problem (25) is equivalent 
to the following dual MAC problem: 

K 

E(m) 
Wir] 



(28) 

subject to ^\.x{Qf^) < AiPi + A2P2 

i=l 

where the noise covariance at the BS of the dual MAC is Ai Ai + A2A2, and r|™'' denotes the achievable 
rate of the iih. user. The problem (28) is a convex optimization problem that can be solved via a standard 
interior point algorithm. With the optimal solution of the problem (28), the optimal solution for the 
problem (25) can be obtained by the MAC-BC transmit covariance matrix transformation. We next 
consider the minimization problem 

min (7(Ai,A2). (29) 

Ai>0,A2>0 

Since the function 5'(Ai,A2) is not necessarily differentiable, we can solve the problem (29) via the 
subgradient algorithm or ellipsoid algorithm. The subgradient of the function 5'(Ai,A2) can be found 
in the following proposition. Refer to Appendix B for the proof. 
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Proposition 5: The subgradient of the function 5'(Ai,A2) at point [Ai,A2] is [-Pi — tr(QAi),P2 — 
tr(QA2)], where Q is the optimal solution of the inner layer optimization problem (28) with Ai = Ai 
and A2 = A2. 

Remark 3: The Lagrangian function of the problem (24) can be written as 

K 

J2 Win - /ii (tr(QAi) - Pi) - /i2 (tr(QA2) - P2) (30) 

i=l 

while the Lagrangian function of the problem (25) can be written as: 

K 

Y^wm - A(Aitr(QAi) + A2tr(QA2) - AiPi - A2P2). (31) 

i=l 

By observing (30) and (31), we can say that the two Lagrange functions are identical to each other 
if we choose jii = AAi and /i2 = AA2. Thus, the auxiliary variables Ai and A2 can be viewed as the 
Lagrange dual variables. 

Since the function 5'(Ai, A2) is a convex function, the subgradient-based algorithm is guaranteed to 
converge to its optimal solution [26]. According to Proposition 4, when the minimum of 5'(Ai,A2) is 
achieved, the optimal solution of the problem (25) is equal to that of the problem (24). 

In summary, the problem (24) is solved through a two-loop algorithm. By exploiting the general 
BC-MAC duality, the inner loop searches the optimal solution of g{Xi, A2), while the outer loop solves 
the 5f(Ai,A2) minimization problem via a subgradient-based iterative algorithm. The convexity of the 
function g{\i, A2) guarantees that the global optimal solution is achieved. 

C. Relationship to Minimax Duality 

The capacity region computation problem with a per-antenna power constraint considered in [5] is a 
special case of the problem (24). By choosing A to be a diagonal matrix with all its diagonal elements 
being zero except the jth diagonal element being 1, the constraint tr(QA) < P can be viewed as the 
power constraint for the jth antenna of the BS. Different from the method discussed in Section IV-B, 
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the dual MAC problem developed in [5] has a minimax form as follows: 



K 



mm max > "^k^k 



K 



(32) 



subject to ^tr(Qf''^) < tr($) 



tr(Q$) < tr($) 



where $ is a diagonal matrix, 




log 



and the diagonal matrix Q is the noise covariance matrix of the dual MAC. The jth diagonal element of 
# is the power constraint for the jth antenna of the BS. Since the noise covariance is also an unknown 
variable, the existing high-efficiency algorithm for the MAC problem cannot be applied. Instead, a new 
interior point method based algorithm is developed in [5] to solve (32). 

The two constraints in (32) have some redundancy, and can be further simplified via the following 
two observations. 

1) The noise covariance constraint tr(Q$) < tr($) holds with equality when the optimal solution is 
achieved. 

2) Given any positive a, if we replace the constraint in (32) with Xlili ^^(Q!'"^'^^) — tttr($) and 
\x{Q<^) < atr($), the optimal value of the problem (32) does not change. 

The observation 1) can be shown by observing that if the constraint tr(Q^*) < tr($) is satisfied 
with inequality, then the minimization part of the problem (32) does not achieve the optimal solution. 
The observation 2) can be proved through the KKT conditions. If we assume that Q* and Q are the 
optimal transmit covariance matrix and the noise covariance matrix of the problem (32), respectively, 
then it is easy to verify that aQ* and aQ satisfy the KKT conditions of the problem (32) after the 
constraints replacement, and the optimal value of the problem (32) does not change. 

Based on these two observations, we can combine the two constraints in (32) into one constraint 
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^iLi < tr(Q^'), and thus the problem (32) is equivalent to the following problem: 

K 

E(mac) 



mm max 

fc=i 



Q Q >0 



(33) 

K ^ 



subject to ^tr(Qf''^) < tr(Q$). 



1=1 



In our derivation in Section IV-B, we can formulate a similar minimax optimization problem of the 
dual MAC. Combining (25), (28) and (29), we have 



mm max > w,r) ' 



K 

(34) 



subject to ^tr(QS"')) < AiPi + A2P2 



where the noise covariance matrix of the dual MAC is A1A1 + A2A2. In the per-antenna power constraint 
scenario, the problem (34) becomes 

K 

(35) 



subject to ^tr(Q;™^) < tr(5$) 



i=\ 

where S = diag(Ai, ■ ■ ■ , AatJ. The problem (35) is identical to the problem (33) by noting that S = Q 
and Q-™'' = Therefore, the problem (32) and the problem (35) are equivalent to each other. 

Although the general BC-MAC duality in Section III and the minimax duality in [10] have substan- 
tially different formulation, the ways by which they handle the multiple linear constraints are equivalent 
essentially. The general BC-MAC duality based method divides the process into two steps: dual MAC 
problem solving, and multiple constraints handling; while the minimax duality combines the two steps 
together. Essentially, only the dual MAC problem solving step exploits the special BC-MAC reciprocity 
relationship, and the multiple constraints handling step is not specific for the BC problem, i.e., it can 
be applied to solve any optimization problem with multiple constraints. The new method decouples the 
two steps, and thus has more flexibility in exploiting the existing algorithm for a MAC. Moreover, the 
minimax problem (32) is specifically formulated to solve the per-antenna power constraint problem. 
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D. Nonlinear Constraint 

According to the proof of Proposition 1, the general BC-MAC duality requires that the transmit 
covariance matrix of the BC be subject to a linear constraint. In this subsection, we consider the case 
with a convex but nonlinear constraint. The capacity computation under a nonlinear constraint can be 
formulated as follows: 

K 

max y WiVi 

tr (36) 

subject to f{Q) <0 

where the function f{Q) is a general nonlinear convex function. Since the general BC-MAC duality 
does not hold generally under the nonlinear constraint, the problem (36) cannot be solved directly. 
However, as shown in the following proposition, the nonlinear constraint problem (36) is equivalent to 
a single linear constraint problem (refer to Appendix C for the proof). 
Proposition 6: There always exists a linear constraint problem as follows 

K 

max WiTi 

tr (37) 

subject to tr(AQj) < 

with A denoting a constant matrix such that it has the same solution as the original problem (36). 

Remark 4: Proposition 6 illustrates that the nonlinear non-convex optimization problem (36) can be 
transformed into an equivalent linear constraint problem (37), which can be solved by making use of 
the general BC-MAC duality in Section III. However, according to the proof of Proposition 6, the 
parameter A in problem (37) cannot be obtained without the optimal solution Q* of the problem (36). 
Note that it is impossible to obtain Q* before solving this problem, and thus the problem (36) cannot 
be solved by using its equivalence with the problem (37). In the following, we present an example to 
illustrate an iterative algorithm that can find a set of linear constraints, and these linear constraints can 
be used to approximate the original nonlinear constraint. 

Example 1: For simplicity, we consider the MIMO BC with K = 2, Nt = 2, and A^^ = 2. The 
transmit covariance matrix is subject to a nonlinear constraint: (ir{QAi)Y + (tr(QA2))^ < P, where 
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1 

Ai = , and A2 = 

J 

can be formulated as follows: 







1 



. Thus, the capacity computation under the nonlinear constraint 




(38) 



As shown in Fig. 2, by defining pi := tr{QAi) and p2 := tr(QA2), the feasible region R : {pi,P2\pi + 
P2 ^ Pi ^ 0,p2 > 0} is a quarter circle in the nonnegative orthant. According to Proposition 6, the 
problem (38) is equivalent to the following problem 



where q = tr(Q*Aj), i = 1,2, Q* is the optimal solution of the problem (38), and C3 = P. 

To present the iterative process, we use a graphical illustration as shown in Fig. 2. Firstly, we 
arbitrarily select a point 'a' on the boundary of the nonlinear region R, and draw a tangent line to 
the pI + P2 = P curve through the selected point. The tangent line corresponds to a linear constraint 
C*^^^ : c[^''tr(QAi) + C2^^tr(QA2) = cf\ where the superscript '1' denotes the index of the linear 
constraints. The weighted sum rate maximization problem with the constraint C*^^^ can be solved through 
the general duality. Assume that Q*-^-* and r^^^ are the corresponding optimal transmit covariance matrix 
and the obtained optimal weighted sum rate, respectively. Since the feasible region of the problem (38) 
is a subset of the feasible region of C^^\ r^^^ is an upper bound on that of the original problem (38). The 
optimal transmit covariance matrix Q^^\ corresponding to the point b in Fig. 2, is not feasible for the 
original problem. Secondly, we find the point 'c' on the boundary of the region R to be closest to the 
point 'b', and draw a new tangent line through the point 'c'. The new tangent line corresponds to a linear 
constraint C''^^ : c[^''tr(QAi) + C2^^tr((5A2) = cf\ By solving the weighted sum rate maximization 
problem with two constraints C^^^ and C''^^ a new optimal weighted sum rate r^"^^ is obtained, where 
^(1) > ^(2) j]^e iterative process will continue until (tr(Q("^ Ai))^ + (tr((5("^A2))^ < P + e holds. 



2 




(39) 



subject to citr(QAi) + C2tr(QA2) < C3 
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where Q*^"-* denotes the optimal solution of the nth iterative step, and e denotes the prescribed accuracy 
requirement. The obtained r'^") forms a non-increasing sequence with the lower bound on the optimal 
solution of the problem (38). Thus, the algorithm will converge to the optimal solution. 




(IV) Pi 



Fig. 2. The iterative process for the nonlinear constraint problem. (I) Select an initial point 'a', plot a tangent line 'AB', and solve the 
optimization problem with a single linear constraint with respect to 'AB', where point 'b' corresponds to the optimal solution. (II) Find 
a point 'c', which is closest to 'b', on the boundary of the feasible region. (Ill) Plot a tangent line 'CD', and solve the optimization 
problem with the constraints 'AB' and 'CD', where 'd' is assumed to be the optimal solution. (IV) So on and so forth, until the optimal 
solution is achieved. 
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V. Beamforming Problem 

In the preceding section, the capacity computation for the MIMO BC is considered. In this section, 
we consider the MISO BC from a beamforming perspective^ Depending on the objective function and 
the constraints, the beamforming optimization problems can be divided into two classes. One is the 
SINR balancing problem [15], [27], i.e., maximizing the minimum SINR among all the users. The 
other one is the power minimization problem with SINR constraints [5], [24], [28], i.e., minimizing 
some power function with SINR constraints. In the case of Nr = 1, the SINR of the ith user under 
the nonlinear encoding and decoding scheme can be expressed as follows: 

SINRi = ^- ^ ' „ (40) 

where hi denotes the Nt x 1 channel vector from the BS to the ith user, and Ui denotes the transmit 
beamforming vector. Under the linear encoding and decoding scheme, the corresponding SINR can be 
expressed as follows: 



SINRf^ = '^^^^ . (41) 

It has been shown that the SINR constraint can be transformed into semidefinite programming (SDP) 
or SOC form [24], and the corresponding beamforming problem can be solved via the standard interior 
point algorithm. However, the standard algorithm does not exploit the special structure of the problem, 
and may be computationally expensive. A number of efficient iterative algorithms have been proposed 
to solve the beamforming problem. In this section, we will combine these iterative algorithms and 
the general BC-MAC duality to solve the beamforming problem with multiple linear constraints. The 
methods discussed in this section are applicable to both SINR definitions (40) and (41). 

'There are two reasons that we do not consider the MIMO BC scenario. 1) For the MIMO BC case, due to the multiple data streams, 
to impose the SINR constraint for each data stream is only appropriate when the independent encoding for each data stream is applied. 
In the previous work, the beamforming problems are considered only under the MISO scenario [15], [24]. 2) According to the definition 
(7), the SINR of the BC is neither convex nor concave with respect to v and u when A'^t > 1. Therefore, the MIMO optimization 
problem with SINR constraints is not a convex problem, and cannot be solved efficiently. 
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A. SINR Balancing Problem 

In this subsection, we apply the general BC-MAC duality in Section III to solve the SINR balancing 
problem with multiple linear constraints: 

max a 

Q>0,a>0 

subject to SINRj > a^i, V i 

(42) 

tr(QAi) < Pi 
tr(QA2) < P2. 

Note that Q = J2i^=i Pi'^i'^f '^^ this case, and thus the problem can be viewed as a joint beamforming 
and power allocation problem. The SINR balancing problem for the BC has been studied in [15], [27]. 
However, due to the limitations of the conventional BC-MAC duality, previous results can be applied 
only to the case in which there is a single sum power constraint at the BS. Furthermore, the SINR 
balancing problem for the MAC with multiple constraints has been studied in [20], where it is shown 
that the multiple constraints can be completely decoupled. However, the decoupling property does not 
hold for the BC scenario. 

To solve the problem (42), we first consider the following problem: 

fi'baKAi, A2) := max a 

subject to SINRi > 07^, V i (43) 

Aitr(QAi) + A2tr(QA2) < AiPi + A2P2 
where Ai and A2 are auxiliary variables. Similar to Proposition 3, we can prove that the optimal solution 
of the problem (43) is an upper bound on that of the problem (42). The problem (43) can be transformed 
into its dual MAC problem via the general BC-MAC duality and efficiently solved by the iterative 
algorithm in [15]. Moreover, the minimization problem minAi,A2fi'bai(Ai, A2) can be solved through the 
subgradient algorithm or ellipsoid algorithm. The convexity of the function (^ba^Ai, A2) guarantees the 
convergence of the subgradient-based algorithm. Similar to the capacity region computation problem 
in Section IV-B, it can be proved that the algorithm converges to an optimal solution of the problem 
(42). 
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Remark 5: The method to process mukiple linear constraints in the SINR balancing problem is 
identical to that in the capacity computation problem. We first introduce an upper bound function, the 
solution of which can be obtained via the general BC-MAC duality. We next compute the minimum 
value of the upper bound function via a subgradient-based algorithm. Note that the iterative algorithm 
in [15] cannot be applied to the minimax duality case as the iterative algorithm requires the explicit 
expression of the noise covariance matrix. From this perspective, the general BC-MAC duality has 
broader application than the minimax duality. 

B. Power Balancing Problem 

In this subsection, we consider the power balancing problem with SINR constraints. Mathematically, 
the problem is formulated as follows: 

min a 

Q>0,a>0 

subject to SINRj > % 

(44) 

tr(QAi) < aPi 
HQA2) < aP2. 

Since the problem (44) has multiple power constraints, the general BC-MAC duality cannot be applied 
directly. By introducing two auxiliary variables Ai and A2, we transform the problem (44) into the 
following single power constraint problem: 

5'pow(Ai, A2) := min a 

Q>0,a>0 

subject to SINRi > 7^ (45) 

Aitr(QAi) + A2tr(QA2) < a(AiPi + A2P2). 
Similar to Proposition 3, the optimal solution of the problem (45) is a lower bound on that of the 

problem (44). Thanks to the general BC-MAC duality, the MIMO BC problem (45) is equivalent to 

its dual MAC problem as follows: 

K 

min yHQ'r^) 

tr (46) 

subject to SINrJ""^ > 7^ 
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where the noise covariance of this dual MAC is AiAi + X2A2. By solving the problem (46) with the 
algorithm in [28], and utilizing the MAC-BC transmit covariance matrix transformation, the optimal 
solution of the problem (45) can be obtained. Next, we consider a maximization problem as follows: 

max fi-powlAi, A2). (47) 

Ai>0,A2>0 

Similar to Proposition 5, we have the following result concerning the subgradient of the function 

dpowi^l, A2). 

Proposition 7: The subgradient of the function (^powl^i, A2) at [Ai, A2] is [tr(QAi) — Pi,tr(QA2) — 
P2], where Q is the optimal solution of the inner layer problem (45) with Ai = Ai and A2 = A2. 

The proof of Proposition 7 is similar to that of Proposition 5, and thus is omitted here. With 
Proposition 7, the maximization problem (47) can be solved through the subgradient algorithm or 
the ellipsoid algorithm. Similar to Proposition 4, when the maximum of 5'pow(Ai,A2) is achieved, the 
optimal solution of the problem (45) is equal to that of the problem (44) 

Remark 6: The maximum per-antenna power constraint minimization problem was considered in 
[5], which is a special case of the problem (44). In [5], the problem is transformed into its minimax 
dual MAC problem, in which the noise covariance matrix of its dual MAC is an unknown variable. 
A subgradient-based iterative algorithm is developed therein to obtain its optimal noise covariance. 
However, since the noise covariance matrix appears in both the constraints and the objective function, 
it is difficult to have a routine method to obtain its subgradient in [5]. While, in contrast, due to the 
clear physical meaning of the variables Aj (Lagrange dual variables with respect to some constraints), 
the subgradient of the lower bound function (7pow(Ai, A2) can be readily obtained. 

VL Numerical Results 

In this section, we present several numerical results to illustrate the effectiveness of the proposed 
algorithms. For simplicity, we consider a MIMO BC with K = 2, Nt = 2, Nr = 2 for the capacity 
computation problem and with K = 2, Nt = 2, A^^. = 1 for the beamforming problem. The noise 
covariance matrix at each user is assumed to be an identity matrix. 
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A. Capacity Region of the MIMO BC 

In this example, we compute the capacity regions of the MIMO BC with a sum power constraint, a 
per-antenna power constraint and sum power plus per-antenna power constraints, separately. The sum 
power constraint is taken to be 10, and the per-antenna power is taken to be 5. For the sum power plus 
per-antenna power constraint case, the sum power constraint is 8, and the per-antenna power constraint 



is 5. The channel matrices are chosen to be Hi 



with a sum power constraint, the algorithm is simi 



1 
0.2 0.6 



and i?2 



0.5 
0.2 1 



. For the case 



ar to that in [6]. For the case with a per-antenna 



power constraint the subgradient-based iterative algorithm developed in Section IV is applied. For the 
sum power plus per-antenna power constraints case, two different algorithms are adopted. The first 
one is the subgradient-based iterative algorithm. The second algorithm is a heuristic algorithm, and is 
based on the result obtained in the case with a sum power constraint. With the sum power constraint 
solution, the transmit covariance matrix is normalized such that each antenna's power satisfies the 
per-antenna power constraint. The regions obtained by these algorithms are shown in Fig. 3. Since the 
heuristic algorithm obtains the suboptimal solution, the fourth line is just an achievable rate region of 
the MIMO BC with sum power plus per-antenna power constraints. Moreover, since the per-antenna 
power constraint is stricter than the sum power constraint, the capacity region of the case with a sum 
power constraint is larger than that of the case with a per-antenna power constraint. 



B. Weighted Sum Rate Maximization With Nonlinear Constraint 

This subsection is to present the simulation result of the Example 1 in Section IV-D. Suppose that 

0.3 0.2 



2 
0.5 0.6 



and Ho 



1.5 



P in (36) is 100. The channel matrices are chosen to be Hi = 

In Fig. 4 (a), the non-linear constraint function values are plotted versus the iteration steps. It can be 
observed that the non-linear transmit covariance constraint is satisfied when the optimal solution is 
achieved. In Fig. 2 (b), the achieved sum rates are plotted versus the iteration steps. The curve in Fig. 
4 (b) is non-increasing, since the results in the former steps are obtained by solving the weighted sum 
rate problem with relaxed constraints. 
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1.4 




Fig. 3. The capacity regions for the MIMO BC with various power constraints. 



1 
0.5 0.6 



and H2 



0.4 
0.5 1.5 



C. SINR Balancing With Multiple Linear Transmit Covariance Constraints 

In this example, we consider the SINR balancing problem with a per-antenna power constraint. We 
assume that each antenna's transmit power is subject to the constraint 5, and each user's target SINR is 

7, = 1, for i = 1,2. The channel matrix is chosen to be Hi = 

The convergence behavior of the algorithm in Section V-A is shown in Fig. 5. The achieved SINR for 
each iteration is plotted in Fig. 5 (a). It can be observed that the curve in Fig. 5 (a) is non-increasing, 
and the achieved SINR for each iteration is greater than or equal to the final result. This is because 
that the optimal solution of the problem (43) is an upper bound on that of the original problem (42). 
The auxiliary variable values are plotted in Fig. 5 (b). It can be seen from this figure that only one 
constraint is active when the optimal solution is achieved. 
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VII. Conclusions 

In this paper, we have established a general Gaussian BC-MAC duality, where the BC is subject to a 
general linear constraint and the MAC is subject to a weighted sum power constraint. This general BC- 
MAC duality can be applied to solve the capacity computation and beamforming optimization problems 
with multiple convex linear/nonlinear constraints. The relationship between the new method and the 
previous minimax-duality based method has also been discussed. Moreover, it has been shown that, 
compared to the minimax duality the general BC-MAC duality offer greater flexibility for solving BC 
optimization problems. This new duality also generalizes the conventional Gaussian BC-MAC duality. 

Appendix 

A. Proof of Corollary 1: The corollary can be derived from Proposition 1 and the relationship between 
the SINR and the achievable rate (8). We will verify the corollary directly via covariance matrix 
transformation, as follows. 

The achievable rate of the ith user of the dual MAC can be written as 



(m) 1 

r • = log 



log 



log 



(48) 



1 tr-ff^M 



(49) 



(50) 



where the eigenvalue decomposition of A^^ is U^\U, and A^^^^ = A}^"^. According to (50), 
the MAC can be viewed as a virtual MAC with (i/^A"^/^)^ being its channel matrix and its noise 
covariance matrix being an identity matrix. By exploiting the BC-MAC covariance algorithm in [12], 
the achievable region of the virtual MAC is equal to the virtual BC with HkA ' being its channel 
matrix. Thus, the achievable rate of the i user of the virtual BC can be written as 

= log J ^^=' Z^i-^ ^ (51) 

where ^[1, tr(Qfc) = P. By defining = A~^'^Q^,{A~^I^)^ , we have Ef=i tr(AQf ) = P. Thus, 
the achievable rate of the dual MAC with sum power constraint P is also achievable for the primal 
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BC with the constraint tr(AQ) < P. Similarly, we can prove that the achievable rate of the primal BC 
with the constraint tr:(AQ) < P is also achievable for the dual MAC with the sum power constraint 
P. The proof follows. ■ 

B. Proof of Proposition 5: According to the definition of the subgradient, if [si, S2] is the subgradient 

of fi'(Ai,A2) at point [Ai,A2], then we have ^(Ai,A2) > 5'(Ai,A2) + [^1,52] ■ ([Ai,A2] - [Ai,A2])^ for 
any [Ai, A2]. 

The Lagrange function of the problem (25) can be written as 

K 

L{Q, A) = Y,^i^i - A(Aitr(QAi) + A2tr(gA2) - AiA - A2P2). (52) 

i=l 



Thus, the corresponding dual problem is 



min maxLfQ, A, Ai, A2). (53) 

A>0 Q 



We have 



(7(Ai,A2)-(7(Ai,A2) (54) 
=^(Q>A)|J;f? .-L(Q,A)|J;f? (55) 



\Q=Q,X=\ ^V-*' 'V|Q=Q_A=A 

\q=q,x=\ ~ ^y^^ ^)\q=q,x=\ 



>mA)|Jll'A, .-L(Q,A)|J;f? (56) 



= - A(Ai(tr(QAi) - Pi) + AsKQAs) - P2)) + A(Ai(tr(QAi) - Pi) + A2(tr(gA2) - P2)) (57) 

=(tr(QAi) - Pi)(-AAi + AAi) + (tr(QA2) - P2)(-AA2 + AA2) (58) 

=(tr(QAi) - Pi)(-AAi + AAi - AAi + AAi) + (tr(QA2) - P2)(-AA2 + AA2 - AA2 + AA2) (59) 
=(tr(QAi) - Pi)(-AAi + AAi) + (tr(QAi) - Pi)(AAi - AAi) 

+ (tr(Q A2) - P2) (-AA2 + AA2) + (tr(Q A2) - P2) (AA2 - AA2) (60) 

=A[Pi - tr(QAi), P2 - tr(gA2)] ■ [Ai - Ai, A2 - A2]^ (61) 

where (61) is due to (tr(QAi) - Pi)Ai + (tr(QA2) - P2)A2 = 0. The proof follows. ■ 
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C. Proof of Proposition 6: Suppose that Q* is the optimal solution of the problem (36). The KKT 
conditions of the problem (36) can be written as 



dQi 



Q.=Q* 



q.=q: 



(62) 



XfiQl = 0. (63) 



Note that Q = Y.t=iQi^ ^^'^ ^us we have /'(QJ = f'{Q)§Q- = f'{Q)- Now, let us consider the 
linear constraint problem 



K 

max > Wiri 

Q (64) 



subject to tr(AQ) < 
where A = f'{Q) ^ ^ . The KKT conditions for the problem (64) are 



q.=q: ^ ' 



+ (65) 

q.=q: 



dQi 

A/(Q*) = 0. (66) 

It is easy to observe that Q* satisfies the KKT conditions of the problem (64). Combining this with 
Proposition 2, we can conclude that Q* is the optimal solution of problem (64). The proof follows. ■ 
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2 3 

Iteration steps 



(a) The value of (tr(QAi))^ + (tr(QA2))^ versus the iteration steps. 




3.05 
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Iteration steps 



(b) Achieved sum rate versus the iteration steps. 



Fig. 4. The convergence behavior of the subgradient-based algorithm for the weighted sum rate maximization problem with a non-linear 
constraint. 
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(a) Auxiliary variable values versus the iteration steps. 
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(b) Achieved balanced SINR versus the iteration steps. 



Fig. 5. The convergence behavior of the subgradient-based algorithm for the SINR balancing problem with per-antenna power constraint. 
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